Identifying rare variants for genetic risk through a combined pedigree and phenotype approach: application to suicide and asthma.

Suicidal behavior is a complex disorder, with evidence for genetic risk independent of other genetic risk factors including psychiatric disorders. Since 1996, over 3000 DNA samples from Utah suicide decedents have been collected and banked for research use through the Utah Medical Examiner. In addition, over 12,000 Utah suicides were identified through examination of death certificates back to 1904. By linking this data with the Utah Population Database, we have identified multiple extended pedigrees with increased risk for suicide completion. A number of medical conditions co-occur with suicide, including asthma, and this study was undertaken to identify genetic risk common to asthma and suicide. This study tests the hypothesis that a particular comorbid condition may identify a more homogeneous genetic subgroup, facilitating the identification of specific genetic risk factors in that group. From pedigrees at increased risk for suicide, we identified three pedigrees also at significantly increased familial risk for asthma. Five suicide decedents from each of these pedigrees, plus an additional three decedents not from these pedigrees with diagnosed asthma, and 10 decedents with close relatives with asthma were genotyped. Results were compared with 183 publicly available unaffected control exomes from 1000 Genomes and CEPH (Centre d'etude du polymorphisme humain) samples genotyped on the same platform. A further 432 suicide decedents were also genotyped as non-asthma suicide controls. Genotyping was done using the Infinium HumanExome BeadChip. For analysis, we used the pedigree extension of Variant Annotation, Analysis and Search Tool (pVAAST) to calculate the disease burden of each gene. The Phenotype Driven Variant Ontological Re-ranking tool (Phevor) then re-ranked our pVAAST results in context of the phenotype. Using asthma as a seed phenotype, Phevor traversed biomedical ontologies and identified genes with similar biological properties to those known to result in asthma. Our top associated genes included those related to neurodevelopment or neural signaling (brain-derived neurotrophic factor (BDNF), neutral sphingomyelinase 2 (SMPD2), homeobox b2 (HOXB2), neural cell adhesion molecule (NCAM2), heterogeneous nuclear ribonucleoprotein A0 (HNRNPA0)), inflammation (free fatty acid receptor 2 (FFAR2)) and inflammation with additional evidence of neuronal involvement (oxidized low density lipoprotein receptor 1 (OLR1), toll-like receptor 3 (TLR3)). Of particular interest, BDNF has been previously implicated in both psychiatric disorders and asthma. Our results demonstrate the utility of combining pedigree and co-occurring phenotypes to identify rare variants associated with suicide risk in conjunction with specific co-occurring conditions.

Genetic risk factors in two Utah pedigrees at high risk for suicide.

We have used unique population-based data resources to identify 22 high-risk extended pedigrees that show clustering of suicide over twice that expected from demographically adjusted incidence rates. In this initial study of genetic risk factors, we focused on two high-risk pedigrees. In the first of these (pedigree 12), 10/19 (53%) of the related suicides were female, and the average age at death was 30.95. In the second (pedigree 5), 7/51 (14%) of the suicides were female and the average age at death was 36.90. Six decedents in pedigree 12 and nine in pedigree 5 were genotyped with the Illumina HumanExome BeadChip. Genotypes were analyzed using the Variant Annotation, Analysis, and Search program package that computes likelihoods of risk variants using the functional impact of the DNA variation, aggregative scoring of multiple variants across each gene and pedigree structure. We prioritized variants that were: (1) shared across pedigree members, (2) rare in other Utah suicides not related to these pedigrees, (3) < or = 5% in genotyping data from 398 other Utah population controls and (4) < or = 5% frequency in publicly available sequence data from 1358 controls and/or in dbSNP. Results included several membrane protein genes (ANO5, and TMEM141 for pedigree 12 and FAM38A and HRCT1 for pedigree 5). Other genes with known neuronal involvement and/or previous associations with psychiatric conditions were also identified, including NFKB1, CASP9, PLXNB1 and PDE11A in pedigree 12, and THOC1, and AUTS2 in pedigree 5. Although the study is limited to variants included on the HumanExome BeadChip, these findings warrant further exploration, and demonstrate the utility of this high-risk pedigree resource to identify potential genes or gene pathways for future development of targeted interventions.